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THE PALM MEASURE AND THE VORONOI TESSELLATION 
FOR THE GINIBRE PROCESS 

By Andre Goldman 

Universite Claude Bernard Lyon 1 

We prove that the Palm measure of the Ginibre process is ob- 
tained by removing a Gaussian distributed point from the process 
and adding the origin. We obtain also precise formulas describing 
the law of the typical cell of Ginibre- Voronoi tessellation. We show 
that near the germs of the cells a more important part of the area 
is captured in the Ginibre- Voronoi tessellation than in the Poisson- 
Voronoi tessellation. Moment areas of corresponding subdomains of 
the cells are explicitly evaluated. 

1. Introduction and statement of the main results. The Poisson- Voronoi 
tessellation is a very popular model of stochastic geometry. This is mainly 
due to its large range of applicability: crystallography [7], astrophysics [32] 
and telecommunications [1], to mention only a few. This is also due to the 
simplicity of the simulation procedures [11, 13, 31], and to the fact that sev- 
eral theoretical results related to its geometrical characteristics are available 
[4, 9, 10, 21, 22]. An extensive list of the areas in which these tessellations 
have been used can be found in [23, 29]. Nevertheless, the other side of the 
picture is that the comparative triviality of this model makes it inappro- 
priate to describe precisely some natural phenomena. Hence, it seems both 
interesting and useful to explore other random point processes and their 
Voronoi tessellations. For instance, Le Caer and Ho [16] describe, by means 
of Monte Carlo simulations, statistical properties of the Voronoi tessella- 
tion associated to the Ginibre process of eigenvalues of random complex 
Gaussian matrices [8] (see also [26]). The idea behind their study is that 
the repulsive character of the distribution of random points makes the cells 
more regular. Consequently, the associated tessellation fits better than the 
Poisson- Voronoi one, as in, for example, the structure of the cells of biolog- 
ical tissues. We recall that the Ginibre process [20, 28], is a determinantal 
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process <f> C R , both isotropic and ergodic with respect to the translations 
of the plane R 2 = C, with the integral kernel, 

(1) ^ 1 ,z 2 ) = (l/7T)e^ 2 exp(-(l/2)(|z 1 | 2 + |z 2 | 2 )), ( Zi ,z 2 )gC 2 . 

It is also pertinent to consider the full class of determinantal processes 4>* a re- 
lated to the kernels, K* a ( Zl ,z 2 ) = (l/vr)^ 1 /"^ 2 exp(-(l/2a)(|zi | 2 + \z 2 \ 2 )), 
(z 1 ,z 2 ) G C 2 , with < a < 1. The process (f)* a can be obtained by deleting, 
independently and with probability 1 — a, each point of the Ginibre process 
(j) and then applying the homothety of ratio \fa to the remaining points in 
order to restore the intensity of the process (p. Besides, it is easy to verify 
that 4>* a converges in law when a — > to the Poisson process. In other words, 
the processes 4>* a constitute an intermediate class between a Poisson process 
and a Ginibre process. In order to challenge the classical Poisson- Voronoi 
model, it is necessary to have some theoretical knowledge about geometric 
characteristics of Ginibre- Voronoi tessellations. The main tool for this is the 
notion of a typical cell in the Palm sense [22]. To explain this notion, we 
introduce, for a general stationary process ip, the notation, 

^o = (^|OG^)\{0}. 

The typical cell of ip is 

C = {z G C; Vu G ipo, \z\ < \z — u\}. 

When ip is ergodic, the laws of the geometric characteristics of the typical 
cell coincide (see [3, 5, 6]) with the empirical distributions of the correspond- 
ing characteristics associated to the Voronoi tessellation, {C(u,ip);u G tp}, 
whose cells are 

C(u,tp) = {z G C; V-u G i/j, \z — u\ <\z — v\}, u Eift. 

If -0 is a Poisson stationary process, then the Slivnyak formula [22] states 
that, for every finite set ScC, 

(4>\Sci>)\s l = 4>. 

Hence, in this case, the Palm measure of ip is the law of the process ip U {0} 
obtained by adding the origin to ip. For every determinantal process tp, a 
result obtained by Shirai and Takahashi [27] states that ipo is determinantal 
as well. It follows that in the Ginibre case the process 4>o is determinantal 
with the kernel, 

K ( Zl ,z 2 ) = (l/7r)(e*« - l)exp(-(l/2)(|z 1 | 2 + |z 2 | 2 )), 

(2) 

{ Zl ,z 2 )eC 2 , 

and that 

C l = C(O,0 O ). 
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Note that the process cfto is nonstationary. 

Our first main result is that 4>o can be obtained from <j>, simply by deleting 
one point. A more precise statement follows. 

Theorem 1. There exists a Gaussian-centered random variable Z , such 
that E\Z\ 2 = 1 and 

0^0oU{Z}, fon{z} = 0. 

Theorem 1 tells us that there exists a version of the Ginibre process <fi 
such that the Palm measure of is the law of the process obtained by 
removing from <p a Gaussian-distributed point and then adding the origin. 
As an intermediate step on our way to further results, consider a locally 
compact Hausdorff space E with a countable basis and a reference Radon 
measure A, and a general stationary determinantal process !/>cE with kernel 
K defined on E 2 . We introduce the following conditions. 

Condition I. The measure A has full support, that is, for every open 
set U C E, X(U) is positive. 

Condition A. The kernel K is a continuous function on E 2 . 

Condition B. For every bounded Borel set A C E, all the eigenvalues 
of the operator Ka, acting on L 2 (A, A), lie in the interval [0, 1[. 

Our intermediate result is the following. 

Theorem 2. Assume that E and i/j satisfy Conditions I, A and B. Then 
the process ip® is stochastically dominated by the process ifj. 

More generally we prove that: 

Theorem 3. Assume that E satisfies Condition I and consider two ker- 
nels K and L satisfying Conditions A and B above. Denote by tp the deter- 
minantal process associated to the kernel K and by ip the process associated 
to the kernel L. Suppose that K > L in the Loewner order. Then the process 
ip dominates stochastically the process ip. 

Recall that K > L in the Loewner order if K — L is a positive semidefinite 
operator. For the kernels K and Kq defined by formulas (1) and (2), we have, 
obviously, K > Kq thus Theorem 3 implies Theorem 2. 

The proof of Theorem 3 rests on an explicit description of the marginal 
laws of the process ip, obtained in Section 1, which allows us to use a similar 
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result proved by Lyons [19] (for commuting operators) and for Borcea, Bran- 
den and Liggett [2] (without this restriction) in the discrete determinantal 
process setting (see also errata to [19] on Russel Lyon's website). 

Thanks to the characterization of the Palm measure, we obtain, following 
[4], precise formulas (see Section 5) which describe the law of the typical 
cell of the Ginibre-Voronoi tessellation. The integrals involved are rather 
awkward; this should not be a surprise, since this is already the case for the 
Poisson-Voronoi typical cell [4]. 

In the last part of this work we compare the moments of the areas of 
Poisson-Voronoi and Ginibre-Voronoi cells. We show that near the germs 
of the cells, a more important part of the area is captured in the Ginibre 
case; farther from the germs of the cells, the situation is reversed. That is, 
roughly speaking, Ginibre cells are more stocky than Poisson cells. 

To be more precise, we introduce some notation. Let C p denote the typical 
cell of the Voronoi tessellation associated to a stationary Poisson process in 
C with the same intensity as the process (j). For every positive r and every 
z G C, let B(z,r) C C denote the disc centered at z with radius r, and 
B(r) = B(0,r). For every finite set S C C, 

D(S)=[jB(z,\z\). 

Let V(S) denote the area of D(S). For every Borel set A C C, and every 
positive integer k introduce 



V k (A) 



k 

dz 

2 



where, for z = x + \y in C with (x, y) E R , one sets dz = dx dy. Finally, for 
z = (zi, ... , Zk) € C fc , one sets dz = dz\ • • • dz^. 

Theorem 4. Let k denote a positive integer. 
(a) When r — > 0, 

(3) EV k (C n B{r)) = EV k (C p n B(r))(l + r 2 W k + o(r 2 )) 

with 



W h = -ml V(z)dz. 



7T J B {l)h 

(b) For every positive R, 
(4) EV k {C \ B(R)) < EV k (C p \ B(R)) • e ( 3 / 2 )~ J ( R ) 

with 



J(R) = T L [ e-^- z ^dz 1 dz 2 . 
27r Jb(r)z 
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Hence, there exists a positive constant c such that, for every positive R, 
EV k (C \ B(R)) < EV k {C p \ B(R)) • e (3/2) " cfl2 . 

We are also interested in the location of the point Z in Theorem 1 with 
respect to the process fa. To describe this, let M(4>) denote the set of points 
z G (j) such that the bisecting line of the segment [0, z] intersects the boundary 
of the cell C(0,(p) where we recall that 

C(0, 4>) = {z£ C; Vn £ <j), \z\ <\z- u\}. 

For every set 5cC, let 

(5) H (S) = (7rR D(s) (0,0) - 1) H(l - a n (S)), 

n>0 

where Rd(s) ^ s the resolvent kernel, and a n (S) for n > 0, are the eigenvalues 
of the integral operator K acting on the space L 2 (D(S),dz). 

Theorem 5. For every positive integer k, 

H(z) dz 

(6) P{Z£AT(<P)}> 



f ck yjH{z) dz 



Taking k = 1, Theorem 5 yields 



(7) 

l-/ +00 n„>o( r (n+l,t)/n!)^ 



> 



Jo + °° V^n>i(* n e-Vr(n + 1, *))) n„>o( r (^ + 1, *)/«') dt 
where T(n,t) denotes the incomplete gamma function, defined as 

/+oo 
e~ u u n ~ x du. 

This implies the simpler bound, 

P {ZeN{4>)}>± 

Section 2 contains the necessary background; the key results are Propo- 
sitions 3 and 12. In Section 3 we prove Theorem 3. Theorem 1 is proved in 
Section 4 as a consequence of Theorem 2, Strassen's classical result and be- 
cause of the fact that the radial processes \cj>\ and |0o| are explicitly known. 
Unfortunately, the correlation between the process and the random point 
Z is still unknown. Nevertheless, Theorem 5 gives some partial insight in 
this direction. Finally, we mention that we state our results for the Ginibre 
process, but that it is easy to deduce the corresponding formulations for the 
processes (j)* a with < a < 1 . 
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2. Preliminaries. Let tp denote a point process [28] on a locally com- 
pact Hausdorff space E with a countable basis. For every integer k > 1, let 
ifj( k ) — {z C ip; \z\ = k}, tp^ = tp, be the associated fc-dimensional process. 
Let Hk denote the corresponding intensity measure. This measure is defined 
as follows. Fix a set z G and consider an arbitrary order z = {z±, . . . , Zk}- 
For every permutation a G pk of the index set {1, . . . , k}, denote 

z a = {z a{1) ,...,z a{k) )eE k . 

For every Borel set A G £>(E fc ) of the space E fc , the sum J2aep k ^a(z u ) does 
not depend on the particular ordering of the set z. By summing for 
and taking the expectation, we obtain 

M{A)=yE e E 1 ^)- 

Consider the space Ai a (E) of the counting measures £ on E such that £(^4) 
is finite for all bounded (relatively compact) Borel sets A C E, and let J- be 
the smallest <r-algebra on _M<j(E) for which the map £ i— > £(A) is measurable 
for every bounded Borel set 4cE. 

The point process tp can be thought of as the random measure 

with values in the measurable space F). Note that the space Ai CT (E) 

endowed with vague topology is a Polish space and that the associated Borel 
o"-algebra coincides with the <r-algebra F (see [17] and [18]). 
For every k > 1, the Campbell measure Ck on E fe x .A/f^E) is 

C fc (M) = l J E7 E E 1 ^'^)' M G £>(E fe ) ® F, 

where, as above, the sum ^2 a&Pk lM(z a ,ip), z a = (z a ^, . . . , z a ^), does not 
depend on the particular ordering z = {zi, . . . , z k } of the set z. 

The disintegration of with respect to the measure /i*. gives, for fj, k 
almost every z = {z\ , . . . , Zk) G E fc , the law of the conditioned process | 5 = 
{z\, . . . ,Z)~} Gip( k ') (see [14]). Campbell formula reads as follows (see [14]). 
Assume that / is a measurable positive function defined on E fc x A^o-(E), 
such that f(z,-) = f(z a ,-) for every permutation a G pk and for fx k almost 
every z G E fc , thus / defines a function acting on sets z G ipw by f(z, •) = 
f((zx, . . . , 2%), •) where 2 = {zi, . . . , z^} is an arbitrary ordering. Then, 

E E f&w 

(8) 

= y Ef((z\, . . .,z k ), (ip I {zi, ...,z k }C ip)) dpt k (zi, ■ ■ -)2fc)- 
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Let A denote a Radon measure on E, that is, a Borel measure such that X(A) 
is finite for every compact set icE. In this paper A will be the Lebesgue 
measure on C = R 2 or the standard counting measure on a finite discrete 
set. The point process ip is determinantal if the following properties hold. 

1. For every k > 1, p k is absolutely continuous with respect to the product 
measure X k on E fc , that is, there exists a density p k such that 

dp k = p k d\ k . 

The density p k is called the correlation function. 

2. There exists a kernel K :K x E — > C which defines a self-adjoint, locally 
trace-class operator, such that, for every z = (z\, . . . ,z k ) in E fc , 

(9) Pk{z) = -rjdet(K(zi, Zj))i<ij< k . 

We use Fredholm notation; hence for every k > 1 and every u = (u\ , . . . , u k ) 
and v = (vi, . . . , v k ) in E fc , 

Kh , ---'l k )=det(K(u i ,v j )) 1 < i , j < k . 

y u\ , . . . , u k j 



Furthermore, 



K | u ) =K I ui,...,u k 



VI \Vl,...,V k 



Assume that tp is determinantal. For every k > 1 and every z S E fe such that 
is positive, let ij) z denote the determinantal process with kernel 

K( u ' z ) 

(10) Kz (u,v) = 1 ^-, (u,v)£E 2 . 

With this notation, for every positive integer k and p and every z € E fc and 



(11) K z 



K( z ' v ) 



Note that, if ^ (!!'") is positive, then K(f) is positive and ip Z)V = (tp z )v 

A result of Shirai and Takahashi [27] (see also [19]) ensures that for \x k 
almost every z = (z%, . . . , z k ) S E fc , 

(12) ip z = (tp | {zi,. . . ,z k } C ^) \ {zi, . . . ,z fc }- 

Recall from [22], that if E is a vector space and ip is a process, stationary 
with respect to the translations of E, then the associated Palm measure Q 
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on (M<r(E), F) is denned by 

where A C E is an arbitrary Borel set such that (ii (A) is positive and finite. 
It follows from (8) (see [27]) that the Palm measure of a determinantal 
stationary process ip with kernel K is the law of the process ipo U {0} where 
■00 = (V' I € \ {0} is determinantal with the kernel, 

(13) A (zi,2:2)- K(0~0) ' z 2 ) EL- . 

If -P{"0 7^ 0} is positive, then K(0, 0) is positive. Applying this to the Ginibre 
kernel, one gets (2). 

Note that if ip is a stationary Poisson process, that is, a point process 
with correlation functions Pk = l/k\ satisfying equality (9) for the degener- 
ate, nonlocally trace-class kernel K{z\,z 2 ) = 5 Zl (z 2 ), then applying formally 
the result by Shirai and Takahashi mentioned above, we obtain Slivnyak's 

formula [22], namely the fact that ip Zl Zk l = tp for every positive k and 
every distinct Zj EC 

For every Borel set A, let Na(iP) denote the number of points of ip in A, 
that is, 

N A (iP) = Y,U(z). 

In the following, we assume that Conditions A and B below hold. 

Condition A. The kernel K(zi,z 2 ) is a continuous function of (z\, z 2 ) E 
E 2 . 

Condition B. For every bounded Borel set A C E, the eigenvalues of 
the operator Ka [acting on L 2 (A)] are in the interval [0, 1[. 

For every bounded Borel set A C E, one sets 

K A 2) (z 1 ,z 2 )= [ K(z 1 ,v)K(v,z 2 )d\(v), (z 1 ,z 2 )eE 2 . 
J A 

(n) 

For every n > 3, K A denotes the iterated kernel of Ka, defined as 
K A n) (z 1 ,z 2 )= [ K(z 1 ,v)K A n - 1 \v,z 2 )d\(v), ( Zi ,z 2 )eE 2 . 

J A 

Conditions A and B above imply that the resolvent kernel, 

Ra(z 1 ,z 2 ) = K(z 1 ,z 2 )+Y / K ( a ) (z 1 ,z 2 ), ( Zi ,z 2 )eE 2 , 

n>2 
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Remark 1 . Note that the resolvent kernel is a continuous function of the 
domain in the sense that if (A n ) n >i is a monotonous sequence of bounded 
Borel sets A n C E such that A n f A (and A is bounded) or A n J. A, then 
RA n (zi, 22) — > Ra(zi, z%) A almost surely. 

It is well known that, for every Borel set A C E, the probability of the 
event {N A (ip) = 0} is a Fredholm determinant, namely, 

(14) P{N A (iP)=0} = det(I-K A ). 

More generally Let n > 1 and (Ai)i<i< n denote n disjoint, bounded Borel 
sets of positive measures X(Ai). Introduce 

n 

A=[jA l . 

i=l 

The Laplace transform of the joint law of random variables N Ai ,i = l,...,n, 
is given by the formula 

(15) Eexp(^-Y^tiN A )j=det(I-Kt,A), U G R+,i = 1, . . . ,n, 

where Ki ^ designates the integral operator K acting on the space L 2 (A,dv) 
with dv(z) = Er=i( 1 - e- u )l Al (z)d\{z). Now, (14) implies 

P{N A (ip) =0} =exp|- I K(z,z)d\(z)-J2- [ K { A n) {z,z)d\{; 

I J A n > 2 n J A 

(16) 
and 

(17) p {NaW = 0} = 1 + J2 { -=^ [ K( V )d\ n (v). 

n>l n ' Ja " ^ ' 



The derivation of formulas (14)-(17) can be found in [27]. 

On the other hand, let us recall (see [25]) C. Platrier's classical formula 
from 1937 (established also by I. Fredholm for k = 1), that is, for every 
positive integer k and every z S the relation, 



(18) 
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From (11), (17) and (18), we deduce that, for every positive k and every 
z £ E fc such that K( z z ) is positive, for every bounded Borel set Ac¥,, 

(19) P{N A {$ t ) = 0} = P{N A (iP) = 0} x 



Remark 2. The kernel Ra — K on E 2 is obviously nonnegative. This 
fact, together with relation (19), implies that, for every positive integer k 
and every z£E*, 

(20) P{N A (^ z )=0}>P{N A (iP) = 0}. 

We now establish some useful results. Let n > 1 and (A;)i<j< n denote n 
disjoint bounded Borel sets of positive measures A(-Aj). Introduce 

n 

A={jAi. 

i=l 

The following proposition gives the joint law of random variables (N^KiKn- 



Proposition 3. Consider n>l nonnegative integers fej such that their 
sum k = k\ + • • • + k n is positive. Introduce 

n 

B = JJ^f', M = {(JV^W)l<i<„ = {h)x<i< n }. 
1=1 

Then 

PROOF. Observe that 
(22) P{M} = E f(~ z ^)i 

where the function / is defined as follows: 

n 

f(z, ^) = l(N A ty \ z) = 0) II C 1*1 = fc < 

where the sums run above the following sets: 
n 

Zi = z VI < i < n, \zi\ = hi- 
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Now apply Campbell's formula (8). We obtain 
P{M} = E f&M 



(23) 



Ef({zi,...,Zk},(ip | {zi,...,z k } ail}))dnk{z\,...,Zk) 
P{N A {{^ | {zi, ..,z t }C^)\K.., z k }) = 0} 



x y n c a ^ ^^i' • • • > 

From property (12) we get 



/n 
P{iV A (^) = 0} Y H^i C ^i) d/*fc(*i, .■■,**) 

(24) 

= rrn ' , , / P[Na(iI>z) = 0}dp k ( Zl ,...,z k ), 
\Ai=ik ] - Jb 

where B = nF=i ' ■ 

The last equality above is obtained by counting partitions, noticing that 
P{A^4(^ 2 ) = 0} depends on the set {z±, . . . ,z k } and that the measure dfj, k (zi, 
. . . , Zk) is permutation invariant, that is, we have 

J l D (zi,..., z k ) dfj, k (zi,. ..,z k ) = J 1d(^(i) , • • • , z a{k) ) dfj, k (zi, ...,z k ) 

for every Borel set DcE' and for every permutation ffGpj,. 
Now, inserting formula (9) in (24) above we get 

It remains to apply formula (19) to obtain Proposition 3. □ 

Remark 4. Formula (23) works for any process ip. In particular if we 
take for tp a Poisson process with intensity measure \i then, by Slivniak's 

formula, ip = (ifi \ {zi, . . . , z k } C ip) \ {z±, . . . , z k }, the fcth-order associated 
intensity measure is dfi k (zi, . . . ,z k ) = (l/k\) dfi(zi) • • -dfi(z k ). Consequently, 
for a Poisson process ip, formula (23) above gives the well-known expression 

Ki<n — (h)l<i<n} 

(25) =—^P{ Na (tP) = 0} / dn k { Zl ,...,z k ) 

U k i ] Jb 
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exp(-n(A)) 



8=1 



Consider now u S E such that u) is positive. The process ip u with ker- 
nel JsT n (zi,z 2 ) = (l/K(u,u))[K(zi,z 2 )K(u,u) - K(zi,u)K(u, z 2 )], (zi,z 2 ) G 
E 2 , fulfills similar Conditions A and B. Indeed, if K(z±,z 2 ) is a continuous 
function of (z\,z 2 ) £ E 2 , then /^(zi, Z2) is a continuous function too. For 
every bounded Borel set A C E, denote by c^m (resp. a Uj A,M) the largest 
eigenvalue of the operator Ka acting on L 2 (A) [resp. of the operator K Uj a 
acting on L 2 (A)\. Notice that the kernel 

K(z 1 ,z 2 ) - K u (z 1 ,z 2 ) = (l/K(u,u))K(z 1 ,u)K(u,z 2 ) 

defines clearly a nonnegative operator, and thus K > K u in the Loewner 
order which implies inequality a U) A,M < ctA,M- Consequently, if Condition B 
is satisfied by K then it is satisfied by K u as well. 

Denote by R u ,a the associated resolvent kernel. Applying the relation (18) 
to the kernel K u , one obtains that, for every zGE*, 

-l) n 



(26) 



+E- 

n>l 



-l) n 



n>l 



■;?.! 



z, v 

Z, V 

V 
V 



d\ n (v) 
d\ n (v) 



On the other hand, 
K 

(27) 



u, z 

It, z 



+ 



n>\ 



til 



K 



A" 



u, z, V 



Ru,A 



d\ n (v) 



(-1)* 



n>l 



711 



K 



A" 



d\ n (v) 



Ra 



u, z 
u, z 



From (10), (17), (26) and (27), we deduce that 



(28) R A 



u, z 
u, z 



P{Na(iP)=0} = K( U , U )R u , a ( ~ )P{NaM = 0}. 



Applying formula (21) to the process tp u and using (28) we obtain the propo- 
sition below. 



Proposition 5. Consider n nonnegative integers ki such that k = k\ + 
1- k n is positive. Introduce the set B and the event M u defined as 

n 

B = Y[Af, M u = {(N Al (ipu))i<i<n = (fe)i<i<»}. 
i=i 
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Then, 



(29) P{M u}= ff^W-?, /W"''W). 

Remark 6. Our notation for vectors of indices are such that equation 

(29) holds more generally, for every positive integer p and every u € E p 
such that K(^j is positive, if only one replaces the factor K(u,u) in the 
denominator by K f"j . 

We now state some simple consequences of Propositions 3 and 5. Denote, 
respectively, by < f3 n < 1 and h n , n > 1, the eigenvalues and the eigenfunc- 
tions of operator Ka- We recall that the eigenfunctions, 

(30) h n (z) = — K(z,v)h n (v)dX(v), \\K\\L?{A,d\) = h 

Pn J A 

are well defined and continuous on E and that 

(31) det(I-K A ) = H(l-p n ). 

n>l 

In what follows, we shall always suppose that the eigenfunctions of operators 
are normalized [as in (30)]. 

Assume now that U C E is a bounded open set and that X(V) is positive 
for every open subset V C U. Let < a n < 1 and f n , n > 1, be the eigenvalues 
and the eigenfunctions of the operator Kjj. A standard result (see, e.g., 
Theorem 2 of [30]) asserts the following. 



Lemma 7. For every (21,22) £U xU, 

K(zi,z 2 ) = y^a n f n {z 1 )f n (z 2 ), Ru(zi,z 2 ) = V" - — — f n {zi) Ufa) , 



and i/ie series are absolutely and uniformly convergent for z\ and z 2 in every 
compact subset of U . 

Remark 8. When the kernel K has the form, 

M 

K(zi,z 2 ) = J ^2/3 n h n (zi)h n (z 2 ), {zi,z 2 ) G E 2 , 

n=l 

with functions h n that are continuous on E and orthonormal on a bounded 
Borel set A, then trivially, 

M p n 

Ra{z\,z 2 ) =y n hnfa)h n fa), (zi,z 2 ) GE 2 . 

n=l 
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Consider the case n = 1 in Propositions 3 and 5. It was proved by Hough 
et al. ([12], Theorem 7) that the random variable Nu{ip) has the distribution 
of a sum of independent Bernoulli (oj) random variables. Explicitly, 

k 

(32) p{ Nu ^)=k}=j2 n (i-«n)n^> 

{ni)in<£(rii)i i=l 

where the sum runs over the indices (nj)i<j<fe such that n\ < ■ ■ ■ < rife. Now, 
by (14) and (31), 

P{N V (V) = 0} = det(J - K V ) = H (1 - a n ) 

n>l 

and thus formula above can be written in the following form: 

k 

(33) P{NuW = k} = P{NuW = 0} £ J] r^-- 

/ s 1 ™i 

(n»)i 1 

Assume that u£U and that i^(w, u) is positive. Using (29) we get 

(34) fww=t}= f| y:\ 

il (u, It) 

where we introduce 

and where each last sum runs over the indices (nj)i<i<fc such that n\ < ■ ■ ■ < 
rik and rii^n for every 1 < % < k. 

Indeed, fix M > 2 and consider the kernels 

M 

Ku,m(Z!,Z 2 ) = y2®nfn{zi)fn{z2) 
n=l 

and 

M 

Ru,M(zi,Z2) = y~]z —fn(zi)f n (z2), (zi,Z 2 )£U X U, 

n=l 

where the functions f n are orthonormal on U. We have 
Ru,Mh Z z )d\ k (z) 

jjk \u,zj 

M k 



(35) =E/»(«)t^-EIIt^- 

ra=l (n 4 ) 4 1 
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e Lnw^ det 



/ fn(u)---f n (Zk) \ 

/™<r(l) ( U )" ' /™ CT (1) ( Z k) 



d\ k (z), 



V /n<T(fc) ( U )"' fn a{k) ( z k) J 

where sums run over the indices {rii)\<i<k such that 1 < n\ < ■ ■ ■ < rif- < M 
and rii^n for every 1 <i <k. 
Observe that 



JJk 



( fn(u)---f n (z k ) \ 
■/W(i) \ U ) ' ' ' \ Z k) 



3=1 



d\ k (z) = f n {u) 



\fn a(k) (u)---fn a(k) (zk)J 

due to the fact that the functions f n are orthonormal on U (and rii^n for 
every 1 <i <k). 

Letting M — > +oo and applying Lemma 7, (29) and (35), we obtain for- 
mula (34). 

Now, by elementary (but somewhat lengthy) computations, which we will 
not detail, we obtain the following proposition. 

Proposition 9. With the assumptions above, 

p{Nu(ip u ) < k} - p{N V ^) <k}= P{ ^ w r 0} s fc , 



K(u, u) 



where 



n>l n (ni)i 1 

and where each last sum runs over the indices (nj)i<j<fc such that n\ < ■ ■ ■ < 
nk and ni^n for every 1 < i < k . 

Proposition 9 implies the result below. 

Corollary 10. Let U be a bounded open set and let < au < 1 denote 
the largest eigenvalue of the operator Kjj . Then for every nonnegative integer 
k, every positive integer p and every u£U p such that K(^j is positive, 

P{Nu(^ u ) <k}<{\- a M )- p P{Nu(iP) < k}. 

Proof. By induction, if u G U, then Proposition 9 and formula (34) 
imply 

P{Nuti>u) <k}- P{N v (ij>) <k}< a M P{NuM = k}. 
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Therefore, 

P{Nu(^ u ) < k} < (1 - a M r l P{Nu(^) < k}. 

Consider now u= (v,w) EU x U p . Recall that ip u = (ip w ) v and that K > 
K w in the Loewner order which implies inequality a w> M < ajf where a w> M 
denote the largest eigenvalue of the operator K Wj a- Thus 

P{Nu(^ u ) <k}<(l- aM)- l P{Nu(^ w ) < k} 

from which we obtain the announced result. □ 

Remark 11. If u G E and z = (^i)i<i<jv G E^, write R\ (^) for the de- 
terminant Ra(u' Z z ) m which one replaces the terms Ra(u,u) and Ra(zi,u) of 
the first column by Ra(u, u) — K(u, u) and Ra^, u) — K(zi,u), respectively. 
Then, 

P{NaM <k} = P{NaW <k} + = ° } / Rf( U AdX\ z ). 

klK{u,u) J A k \u,zj 

One can prove this formula, from Propositions 3 and 5, by induction. 

A further simple consequence of Proposition 3 and Remark 8 is the fol- 
lowing. Consider another locally compact Hausdorff space E' with reference 
measure A', some bounded Borel nonintersecting sets P>i C E' of positive 
measures A'(-Bj) and a point process ip' C E' with kernel 

M 

L B (zi,z 2 ) = S ^a n g n {z l )g n (z 2 ), {z 1 ,z 2 ) G E' x E', 

n=l 

where < a n < 1 and the functions g n for 1 < n < M are defined on E' and 
are orthonormal on 



N 



B = \jB l 

i=i 

Now, let ip C E be a point process with kernel 

M 

K A (zi,Z 2 ) = y^ J Otnfn(zi)f n (z2), (zi,Z 2 ) G E 2 , 

n=l 

where the functions f n for 1 < n < M are continuous on E and are orthonor- 
mal on 

N 



A=\jAi. 
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Assume that the following holds. 

For every 1 < i < N and 1 < n, m < M, 



(36) / fn(z)f m (z)dX(z)= g n (z)g m (z)d\'(z). 

Then the following proposition holds. 

Proposition 12. With the assumptions above, for every (ki)i, 
P{(N Ai ^))i<i< N = (ki)i<i< N } = P{{N B M))i<i<N = (ki)i<i< N }- 

Proof. Let k = k\-\ — • + k^. When k = 0, the result follows from for- 
mulas (14) and (31). Suppose now that k is positive. By Remark 8 we have 

M 

Ra(zi,Z 2 ) = }--. —fn{zi)f n {z 2 ), (zi,Z 2 ) £ A X A, 

n=l 

and 

M 

E a 
r. — —9n{zi)gn{z 2 )-, (zi,z 2 ) G B x B. 

Thus, for (j G pfc we obtain 

k M k 

j=l ni,...,n fc =l j=l 3 

and 

k M k 

II RB{Zj,Z a{j) ) = Yl II i ~ 9n 3 {Zj)gn a . 1{j) fa). 
3=1 m,...,n k =l]=l 3 

Denote C = fj^i ^ and C = ]Ji=i B i* ■ Formula (36) implies that 

/ II f»i dA ( Zl ) " " " dX ( z ^ 

Jc j=i 

Jc 'j=i 



(37) 

k 



Then, expanding the determinants appearing below and using the point (37) 
above, one gets the equality 



(38) 



This and (21) give the result. □ 
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3. Stochastic domination, proof of Theorem 3. In this section, we assume 
that Condition I stated in the Introduction is satisfied. 

Recall that a point process a G .M CT (E) stochastically dominates a point 
process /3 G Ai a (K) if Ef(a) > Ef((5) for every bounded increasing mea- 
surable function / defined on the space (Ai a (E), J 7 ). It is is well known 
[17] that the point process a G M a (K) stochastically dominates the point 
process (3 G A4 a (¥,) if and only if P{a G .4} < P{(3 G .4} for every de- 
creasing event A G T . Consider elementary decreasing events of the form 
{VI < i < M, Na { < ki} G F where M is a positive integer, hi, 1 < i < M, 
are nonnegative integers and A{ C E are disjoint, bounded, Borel sets. 

Denote by C T the collection of sets which are a finite union of such 
elementary decreasing events. The following lemma provides a useful tool in 
order to investigate stochastic domination properties of point processes. 

Lemma 13. The process j3 is stochastically dominated by the process a 
if and only if, for every A^Td, 

P{aeA}<P{(3eA}. 

Remark 14. The proof of Lemma 13 is standard. Similar characteri- 
zations are described, for example, in [17]; however, as pointed out by Yo- 
geshwaran Dhandapani at ENS-DI-TREC (France), this result is not explic- 
itly enunciaded in [17]. For completeness we sketch the proof of it in the 
Appendix. 

We will now prove Theorem 3. Consider two kernels K and L, satisfying 
Conditions A and B stated in the Introduction, such that L < K in the 
Loewner order. Denote by (p the process with kernels L and by i/j the process 
with kernel K. The idea of the proof is the following. By Lemma 13 we need 
to show that for every A G J~d, 

(39) P{iieA} < P{<peA}. 

Fix the set A G JF^. Applying the inclusion-exclusion principle it is easy to 
see that there exist nonintersecting, bounded Borel sets BjCE, 1 < % < N, 
such that P{ip G .4} (resp. P{ip G A}) can be expressed as a finite sum, up 
to the sign, of terms of the form, 

P{Vi G S,A^(t/0 = ki} (resp. P{Vi G S,N B .(ip) = h}), 

where S C {1, . . . , ./V}. Let U be an open bounded set such that 

N 

UD\jBi; 

i=i 
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denote also, Bq = U\ (Ji^i Bj. 

By Lemma 7 we have the spectral decomposition, 



n{Z 2 ), 

n>l 

Ku(zi,z 2 ) = J ^2a n f ri (zi)f n (z 2 ), (zi,z 2 ) G U. 

n>l 

The fact that K > L in the Loewner order reads. 
For all f£L 2 (U), 

2 

(40) Va„ / / n (z)/MdA( Z ) >V/3 n / g n {z)W)d\{z 



n>l 



n>l 



The inequality above implies that for every n > 1, the function g n is of 
the form # n = X^>i a fe/fe G ^ 2 (^)- Denote g^M = YJk=i a kfk and consider 
the nonnegative kernels, 

M 



Ku,M = y^Qn/n(^l)/n 



*2 



n=l 



and 



Lu,M = Pngn,M{zi)g n ,M(z 2 ), {zi,Z 2 ) G £/, 



n=l 



acting on L 2 (U). 

Note that 1 1 || < ||i^c/|| < 1 and H-L^mII < ll-^t/ll < 1 where || • || denotes 
the supremum (operator) norm. Furthermore, if y( M ) is the subspace of 
L 2 (U) spanned by the functions f n with 1 < n < M, then by (40), for each 
function /GV (M) , 



A I 



(41) 



n=l 



9n,M(z)f(z) dX(z) 



u 



2 M 



n=l 
M 



9n(z)f(z)dX(z) 



V 



< 



yZ° n / fn{z)f{z)d\{; 



Denote by 7 n and and h n the eigenvalues and the normalized eigenvectors 
of the operator Ljj^m [acting on L 2 (U)]. The properties above imply that 



(42) 0< 7 n<l and h n = ^ Wf k G V (M) , 1 < n < M. 



k=l 
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At last, 

M M 

Lu,M = y~]Pn9n,M( z l)9n,M( z 2) = ^2lnhn(zi)h n (z 2 ) , (zi,Z 2 ) G U. 

n=l n=l 

Let <f^ C U and V (A/) C U be the processes associated, respectively, to 
the kernels Lu,m and Ku : m- 

Lemma 15. When M — > oo, 

P{Vi G S,N Bi (i>M) = ki] -> P{Vi G 5,iV Bi (V) = h}, 
P{Vi G S, iV Bi {^ M) ) = h] -> P{Vi G 5, iV Bi (p) = 

Proof. The straightforward consequence of (21), (31) and Lemma 7. 

□ 

It follows from Lemmas 13 and 15 that in order to prove that the process ip 
dominates the process (p it suffices to show that, for every M > 1 inequality 
(39) is unchanged if we replace the terms of the form 

P{Vi G S,NbM = ki] (resp. P{Vi G S,N B .(ip) = ki}), 

by the terms 

P{V* G S, N Bi (^ M )) = ha} (resp. P{Vi G S, N B . (<pM) = ki}). 

To obtain this result we use the fact that the stochastic domination occurs 
in the finite discrete determinantal process setting. See Theorem 6.2 and 
Paragraph 8 of [19], errata to [19] on Russel Lyon's website, and [2]. The 
link between our situation and a discrete determinantal process is given by 
the following lemma. 

Lemma 16. Let Bi denote nonintersecting Borel bounded subsets o/E, 
and let 

N 

U=\jB l . 

i=0 

Consider an orthonormal set of functions {l n ,n = 1, . . . , M} C L 2 (U). Let Mi 
denote the dimension of the subspace Vi C L 2 (U) spanned by the functions 
l n l Bi with 1 < n < M. 

Then, there exists orthonormal vectors z n = (z" Q ^ , . . . , z™^ 
for 1 < n < M , such that the following property holds. 

For every < i < N and every 1 < n, m < M , 

N, 



(43) E^)0>«(i)= / ln(z)l m (z)dz. 

3=1 JBi 
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Proof. Since the sequence (l n ) n is orthonormal, property (43) implies 
that the sequence {z n ) n is orthonormal as well. Introduce an orthonormal 
basis {e % j)i<j<N i of the vector space Vi C L 2 (U). Then, 

K 

7=1 

The sequence defined by zj^(j) = A"^ fulfills property (43). □ 

Consider now the vectors z n ,n = 1,...,N, associated with Lemma 15, 
to the eigenvectors f n , n = 1, . . . , M, of the kernel Ku,m, and introduce the 
vectors, 

M N 

v n = (v? 0) ,...,v? N) ) = J2b n k z k eU CAfi > n = l,...,M, 

k=l i=0 

related to functions h n , n = 1, . . . , M, of (42). Notice that for every < i < N 
and every 1 < n, m < M, 

Mi . 

(44) T, v MiK)ti)= h n (z)h m (z)d\(z). 

Consequently, (v n ) n is a set of orthonormal vectors. 
Consider now the sets 

v 

E j = {(i,j);l<i<jV j }, 0<j<N,E=\jE J , 

7=0 

and the discrete kernels defined on E by, 

M 

K M ((h, ji ), (t 2 , 32)) = X] a « z fu ) (ii) z (i 2 ) O2) 
n=l 

and 

M 

((tl , jl) , («2 , 32)) = X Tn Vf ix ) (j'l ) V^j (j 2 ) • 
n=l 

Inequality (41) implies that Lm < -f^M in the Loewner order. Introduce the 
determinantal process xCi? with kernel Lm, and the process C, d E with 
kernel -Km- Proposition 12 and formulas (43), (44) imply that 

P{Vi G S,iY Bl (^ (M) ) = h} = P{Vi G 5,7V Ei (C) = h} 

and 

P{Vi G S,N Bi (<pM) = ki} = P{V* G S,N Ei ( X ) = h}. 
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Consequently if we replace in formula (39), the terms of the form 
P{V* G S, N Bi (if) = h} (resp. P{Vi G S, N Bi ty) = A*}) 
by the terms 

P{\Ji G S, N Bi & {M) ) = ki} (resp. P{V* G 5, (^ (M) ) = fcj}), 
we obtain inequality 

P{CeA'}<P{ x eA'} 

for a suitable decreasing event A' G J-d(E). The above mentioned result of 
[2] and [19] asserts that this inequality is indeed true and thus the proof of 
Theorem 3 is finished. 

We are now in position to apply the celebrated Strassen's theorem. This 
follows from the fact that the space M a (K) of counting measures endowed 
with the vague topology is a Polish space and its associated Borel cr-algebra 
coincides precisely with the cr-algebra T (see [17] and [18]). 

Theorem 6. With the hypothesis of Theorem 3, there exists a point 
process r\ such that 

ip = ipUr), (pHr] = 0. 

Theorem 3 implies Theorem 2. More generally we have: 

Theorem 7. Let ip be a point process satisfying Conditions A and B of 
the Introduction. For all points u such that K(u,u) is positive, the process 
ip dominates stochastically the process ip u . 

Proof. It is obvious that K > K u in the Loewner order. □ 

Problem 1. Prove Theorem 7 directly from (21) and (29). 

4. Palm measure of the Ginibre process, proof of Theorem 1. Recall 
that the Ginibre process <p C M 2 = C is a stationary, isotropic point process 
satisfying Conditions A and B of the Introduction. The reference measure 
A is the area measure of M 2 , and Condition I is trivially satisfied. Moreover, 
for every integer k > 1 and every set of distinct points {z±, . . . , Zk} included 
in C, respectively, in C\{0}, K {^'"' z ^) is positive, respectively, ^o(^'"''* fc ) 
is positive. 

From formula (13) it follows that the process c/>o = (</> | G </>) \ {0} is 
determinantal with the kernel Kq such that 

K (z 1 ,z 2 ) = (l/7r)(e^ 2 - l)exp(-(l/2)(|^| 2 + \z 2 \ 2 )), (z u z 2 ) G C 2 . 
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The intensity measure /io,i is absolutely continuous with respect to the mea- 
sure A and has for density the correlation function 



(45) K (z,z) = (lM(l-e-\*\ ). 

In particular, the process <f)Q is not stationary. 

Remark 17. The stationarity of the Ginibre process </> is expressed by 
the fact that for each fixed a G C the determinantal point process with kernel 
K such that K(z\, z 2 ) = K(z\ — a, z 2 — a), that is, 



K( Zl , z 2 ) = (l/7r)e (2l - a)(z2 - a) - (1/2)(|zi - a| +|22 " a| >, ( Zl ,z 2 ) G C 2 , 
coincides (in law) with cp. Note that K ^ K. 

Consider now the radial processes \4>\ and \(j>o\- The result below is well 
known [12, 15]. 

Theorem 8 (Kostlan). Let X n ^ m with n > 1 and m > 1 denote i.i.d. 
random variables with exponential distribution e~ x dx on x > 0. For every 
n > 1, let 



Rr, 



\ m=l 



Then, the collection of moduli of the points of (j) has the same distribution 
as the collection of random variables {R n ,n > 1}. 

(46) H= W {#n,n>l}. 

Remark 18. Note that Theorem 8 implies that, almost surely, there 
exists no (zi,z 2 ) G <fi x 4> such that z\ ^ z 2 and \z\\ = \z 2 \. 

We will shown that result (46) can be deduced from formula (15). Indeed, 
let us fix < r\ < ■ ■ ■ < r n = r and consider the sets A\ = B(r±), Ai = {z G 
C; ri-\ < \z\ < rj} for i = 2, . . . , n, B(r) = UILi ^i- Also, let us fix t{ > 0, 
i = 1, . . . ,n. 

Observe that the functions, 

|2 



f n {z) = {l/V^.)e~^ 2 ^ z n , zeB(r),n>l, 
are orthogonal on B(r) with respect to the measure, 



dv(z) = Y,a-e- U )l Al (z)dHz). 



i=l 
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Denote a n = f B r r \ \f n \ 2 dv{z) then normalizing, we obtain 



J(r) 

|2 , u„|2<, 



K( Zl ,z 2 ) = (l/<K)e z ^-W^ +lz2n = ^2a n f n ( Zl )f n (z 2 ) 



n>l 



with f n (z) = (l/y/a^)f n . Consider now the radial process \cj)\ and the inter- 
vals h = [0,ri], ...,/„=] 

r n— ii r n]- Formulas (15) and (31) imply that 



(47) 

= det(J-^) = Jj(l-on). 

n>l 

Computing (what is an elementary exercise) the Laplace transform 
Fexp(— Yl7=i tiNit (7£)) for the point process 1Z = {R n , n > 1} gives exactly 

the same value. Thus \<p\ l = {R n ,n > 1}. More generally, if ip(F), Fcff, 
F^0, is the point process related to the kernel 



K(F)(z 1 ,Z 2 ) = ^ a nfn(zi)f n (z 2 ), 

then 

(48) |0(F)| l ^{R n ,neF}. 

In particular, 

\M l ={Rn,n>2} 

and 

{R n ,n>l} = {R n ,n>2}U{R l } 

provide a disjoint coupling of \4>o\ and {Ri} with union marginal \(f>\. 
Consider now, for M > 1, the kernels 

1. F M (^i^2) = (lA)E^ ((^2r/^)exp(-(l/2)(|zi| 2 + |z 2 | 2 )), 

2. K Mzi,z 2 ) = (lMZnU(^2)Vnl)eMHm(\zi\ 2 + \z2\ 2 )), 

and denote by the point process associated with the kernel Km and 
by 4>q the point process associated with the kernel ifo,M- 
Observe that on one hand we have 

£card{^ (M) }= f K M (z,z)dz = M, 
Jc 

E card{4 A/) } = f K 0tM (z, z) dz = M - 1, 
Jc 
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and on the other hand, by (9), the correlation function of order M + 1 
(resp. M) for the process <fi( M ^ (resp. 4>q^^) is equal to zero which implies 
that caxd{(j)( M }} < M and card{<^, } < M — 1, almost surely. Therefore, 
cardj^^} = M and card{(/>Q A ^} =M — 1, almost surely. Moreover, we have 
Kq,M < Km in the Loewner order. Formula (48) implies also 

(49) \^\ ^ {R n , 1 < n< M + 1}, \</tf\ ^ {R n , 2 < n < M + 1}. 

It follows from the properties above and from Theorem 6 that there exists 
a disjoint coupling, 

(50) ^ = 4 M) U V W , ^ n J*) = , 

such that the point process rf M ' is a single random variable, r/ M ) = {Zm}- 
By equation (49) and the fact that Remark 18 also applies to the process 

we deduce that \Z M \ l = Ri- 

Lemma 19. The random variable Zm is centered Gaussian, and 

E\Z M \ 2 = l. 

Proof. The random variable \Zm | 2 has exponential distribution e~ x dx, 
x > 0, thus it suffices to show that the law of Zm is invariant by rotations 
O with its center at the origin, that is, P{Zm G ^4} = P{Zm G 0(A)}, for 
every such O. Simple computation gives 

(51) P{Z M eA}= [P{NA(4 M) )< k }- p {N A (^ M) )<k}}. 

0<k<M+l 

The processes <p^ M ^ and 0q are isotropic, hence formula (51) implies the 
result. 

Consider now the laws P^ M \ M > 1, of random elements (Zm,^') with 
values in the product space C x Ai a (C) [the space Ai a (C) being endowed 
with the vague topology]. 

Denote, respectively, by Q, Q^ M \ Qo and Qq , the laws of the processes 

</>, (f> {M) , <p and 4 M) . Finally, let J:C x M a (C) — >• M a {C) be the contin- 
uous application defined by I(x, £) = {x} U (. □ 

Lemma 20. The following properties hold: 

1. The sequences (Q( m ^)m an d (Qq )m o,re tight. 

2. QW A Q andQi M) A Q . 

M-H-oo M-t+oo 

3. The sequence (P^ m ^)m is tight. 
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4. Consider the probability p( M ) on (C x M a (C), B(C) ® J 7 ), tfien I = Q (A/) . 

Proof. Property 1 is obvious from the characterization of tightness for 
random measures (see [14], page 33). Property 2 follows from property 1 and 
Proposition 3. Property 3 is a consequence of property 1 and the fact that, 
by Lemma 19, the standard normal law coincides with the marginal law on 
C of the probability p( M ) . Finally, property 4 is nothing but the equality in 
law (50). □ 

It is well known that a suitable subsequence of (P^ M ')m converges in 
distribution to a probability P* on C x M a (C). Lemma 20 implies that P* 
has, for marginal laws, the standard normal law and Qq and that with P* 

on C x _A/fcr(C), we obtain I = Q. Consequently, a random element with dis- 
tribution P* provides a disjoint coupling (Z, fa) of <j). The proof of Theorem 
1 is then finished. 

One can notice also that we have 

(52) P{Z eA} = ]T[P{ A^o) < k} - P{N A (^) < k}]; 

k>0 

thus, if U is an open set containing the origin, then by inserting the formula 
of Proposition 9 in (52), we obtain 

n>l L k>l (m)i 1 

where the last sum is over the integers (nj)i<j<fc such that rtj < n^i and 
rii^n for every i. Hence, 

v ' ' n>l i>l 1 

and, finally, 

1 ' A' (0,0) nj v 
Thus we find again that the law of Z is Gaussian. Notice also the formula, 

(5 3) P{Z€AlNAiM=Q]=1 -im., 

which follows from (19) via the identities, 

P{Z e A,N A {fo) = 0} = P{N A (M = 0} - P{N A {<t>) = 0} 
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and 

p ^«=°>=S p{ ^ (W=o} ' 



Problem 2. This is an open problem, that is, to know how the random 
variable Z is correlated with the point process 4>q. A similar unsolved prob- 
lem arises in the framework of finite discrete determinantal processes (see 
question (10.1) in [19]). 

Remark 21. The method we used to prove Theorem 1, that is, a cou- 
pling result [formula (50)] in a finite-dimensional case associated with a 
"tightness argument" (Lemma 20) is very similar to that used by R. Lyons, 
in the discrete determinantal process setting, to prove Proposition 10.3 in 
[19]. 

Remark 22. Theorem 6 can be applied to the processes cj> and (fio and 
thus provides a disjoint coupling (f> = ipoL) {rj}. However, there is a difficulty 
to deduce Theorem 1 directly from this (due to the fact that it is unclear 
that the process ij could be taken as being a single random variable). 

Remark 23. The random variables R^, n > 1, are Gamma(n, 1) dis- 
tributed and independent. They are stochastically increasing but not almost 
surely increasing. It is interesting to note that if R n = y/Ylm=i n — 1 > 

is the radial process of a Poisson stationary process which has the same in- 
tensity (1/vr) dz as the process (j>, then the random variables R^, n>\, are 
Gamma(n, 1) distributed as well; they are almost surely increasing and (of 
course) is not independent. 

Problem 3. Construct explicitly random variables Z n , n > 1, such 
that: 

1. 4> l = {Z n ,n> 1}; 

2. Vn> l,\Z n \ = R n ; 

3. 4>q 1 = {Z n ,n>2}. 

Remark 24. The Palm measure of <j)* a is obtained by adding the ori- 
gin and deleting the point \faZ if the latest belongs (which occurs with 
probability a) to the process 4>* a . 



Remark 25. Similar results could be proved for the point process in 
the unit disk of C related to the Bergman kernel and studied in [24]. 
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5. Ginibre Voronoi tessellation, proof of Theorems 3 and 4. Consider 
now the space K, of compact convex sets of R 2 = C endowed with the usual 
Hausdorff metric. For every point process tp, let 

C(u, ip) = {z G C; Vf G ip, \ z — u\ < \z — v\}, 

and let {C(u,<j));u G 0} denote the Voronoi tessellation generated by the 
Ginibre process <f>. Recall that its statistical properties, namely its empirical 
distributions (the process being ergodic), are described [3, 5, 6] by the typical 
cell C defined by means of the identity, 

Eh(C) = ^-E Y, HC(z)-z), 

where h runs through the space of positive measurable functions on /C, and 
B C C is an arbitrary Borel set with the finite positive area A(-B). Consider 
now the cell 

C(O,0 O ) = {z G C; Vu G 0o, |*| <\z- u\}. 
Campbell's formula (8) gives the identity, 

Eh{C) = ^-E Y h(C(0^-z))=Eh(C(0,<^)). 

Hence, 

(54) C^ W C(O,0 O ). 

In what follows, we shall use the notation C(O,0o) = C(0)- The law of the 
random set C(0) can be obtained by means of the method described in [4]. 
Let us introduce some notation. Fix k > 1. 

• For every u G C, let H(u) ={z£C;(z- u,u) < 0}. 

• For every z G C k with z = (z%, . . . , Zk), let Sj(z) denote the intersection of 
half-spaces, 

A' 

fi(z) = f)H(zi/2). 

i=l 

• For every z G C fc , let F(z) = U u efl(z) B{u, \u\) where B(u,r) denotes the 
disk centered at u and of radius r > 0. 

• Let A<zC k denote the set of z G C k such that Sj(z) is a bounded polygon 
with k sides. 

Theorem 9. For every k>3, 

P{C has k sides] = y J P{N T(g) (<f>o) = 0}R ^ (z) ( Z \ dz. 
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Proof. Observe that 

P{C(0) has k sides} = E £ x l {0} (JV>( a ,... A )(^o \ *))■ 



With Campbell's formula (8) applied to the process 0o, one gets 

P{C has ft sides} = P{C(0) has ft sides} = / P{N T{z) ((j) 0:Z ) = 0}dn , k (z) 

J A 

and hence by formulas, (9) and (19), we obtain 

P{C(0) has k sides} = I jf P{JN^ W (<fo, z ) = 0}^ f d* 

= ^ jf P{JV> W (<M = 0}iV(z) ( * ) dz - □ 

In the same way, one can compute, conditionally on the fact that the cell 
has k sides, the expectation of an arbitrary, measurable, positive functional 
of C which is expressed through a function / acting on points {z\, . . . , z k } = 
M(4>o) C 4>o f° r which the bisecting line of the interval [0, Zj\ intersects the 
cell C(0). The resulting integral will have the form 



dz. 



(55) ^ J f{z)P{N Hz) (0q) = 0}^*) ( Z z \ 

Deconditioning, one can obtain analytical formulas of the laws of the geo- 
metric characteristics of the typical cell C. Note that formula (17) gives an 
analytical expression of the probability P{Nj:r z \(<f>o) = 0} which appears in 

(55) . Unfortunately, these integrals are complicated and numerical compu- 
tations are difficult. This drawback appears already in [4] for the typical cell 
of the Poisson-Voronoi tessellation. 

A general result asserts [22] that the first-order moment of the area V(C) 
of the cell C is equal to EV(C) = vr. The moments EV k (C) of higher orders 
can be expressed in terms of integrals more tractable than (55). Recall our 
notation 

k 

D(z 1 ,...,z k ) = \jB(z i ,\z i \)cC. 
i=i 

We use the fact that 

(56) zeC(0) k N D(z) (<p ) = 0. 
Let A C C be a Borel set. From (54), (56) and (19), 



(57) 



E[V k (CnA)}= [ P{N D{z) ( ( j )o ) = 0}dz 

JA k 



J Ak if (0,0) P {N D{z) tt) = 0}dz 
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and hence by (16), we obtain 

E[V k (CDA)}= [ expj- f K (u,u)du 



A k I JD{z) 



(58) 



n>2 

and 



~H\J k omz)( u > u) du } dz 



(59) 



n>2 JD ( Z ) > 



where V(z) denote the area of the set D(z) = D(z±, . . . , Zk)- 

We will now use formulas (58) and (59) in order to compare the area of C 
with the area of the typical cell C p of the Voronoi tessellation associated with 
a stationary Poisson process which has the same intensity measure (l/n)dz 
as the process <f>. For every Borel set A, 

(60) E[V k {C p C\A)] = [ e~ v ^^dz. 

JA k 

Hence by (58), (45) and the fact that for every z = (z\,...,Z)-) € B{r) k , 
D(z) = D(z±, . . . , Zk) C B(2r), one has 

(61) E[V k (CDB(r))} < E[V k (C p D B(r))]ex V ( - [ dz) . 

V 71 " JB(2r) ) 

5.1. Proof of Theorem 4, part (a). Formula (60) implies that 

(62) E[V k {C p C\B{r))]= f e~ v ^ ^ dz = r 2k f e -** v W* dz. 

JB{r) k JB(l) k 

For z £ C k , let ao,n, n > 1, denote the eigenvalues of Kq acting on D(z) 
where uq i is the largest eigenvalue, then 



(63) 



m>2 \ ) m >2 n>l 



n) 



2(l-a ,i)^ 



Furthermore, 



^(a ,n) 2 = / ^ 2 i (z) (n,n)dn=(l/7T 2 ) I |1 - e ««|2 e -l«l a -H a dttdt7 . 

n >l ' JD(2) 2 
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If z € B(r) k , then D(z) C B(2r), and hence 

V(a n) 2 <(l/vr 2 ) / \l-e uV \ 2 e^ u \ 2 ^ 2 dudv = 0(r 8 ). 

This and (63) imply (since «o,i decrease when domain decrease) 
(64) Y,^L ) <%)M du = °^)- 



m>2 

Therefore, by (58) we obtain 

E[V k (CnB(r))}= [ expj- f K (u,u)du 

JB(r) k I JD(z) 

(65) -E-/ {u,u)du\dz 

= (l + 0(r 8 )) f exp(- / Ko(ti,u)d«ldz. 

Moreover, 



B(r) k I JD(«) 



/ exp< — / Kq(u, u) du > dz 

JB(r) k [ JD(z) J 

fc / expl- K (u,u)du\dz 



'B(l) fe I Jr_D(z) 

(66) 

„2fc / „, I „2 



r 



/" exp(-r 2 iy(z) + — / e" |rn|2 du) dz 

.7.8(1)* I 71 71 J 



2 

l + -y(z) + 0(r 4 ) 



7T 



exp<j —r 2 —V(z) ^ dz. 



Formula (3) is a straightforward consequence of (62), (65), (66) and the 
asymptotic equality, 

. . f B(1)k V(z)eM-r 2 (lMV(z)}dz x r 

67 r J — 27T7 u// \ ~i 7 = — / Fzdz + Or 2 . 

Jb(i)* exp{-H(l/7r)F(z)}dz vr fc 

5.2. Proof of Theorem 4, part (b). For z G {C\B(R)) k , let a n , /„, n > 1, 
denote the eigenvalues and the eigenfunctions of the operator K on D{z) 
where ot\ is the largest eigenvalue. 

By Lemma 7, 

R D{Z) (0,0) E™>iK/(l-an))|/n(0)| 2 1 

y b °) rFFn~n\ _ ~ I t /n\l2 - 



^(0,0) En>l«n|/n(0)| 2 ~ 1 - <*{ 
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One has 
(69) 
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^2a n = (l/7r)V(z). 

n>l 



Note also that when z € (C\ B(R)) k , there exists a £ C, \a\ = R such that 
D(z)D B(a,R). Thus 



(70) 



y](an) 2 = / K {2) ( ) (u,u)du = \ [ 

^2 J D(z) D[Z) ^Jd(z) 2 



-\ u ~ v \ 2 dudv 



[ e -l«-^l 2 dudv = \ I 

JB(a,R) 2 11 Jb(R) 2 



1 

" Jb(cl,R) 2 



dudv. 



Introduce 



#(0,0) 

Thus, by (14), (31), (68) and (69), 
(*) < expj ^log(l - a n ) \ < exp 

and hence by (70), we obtain 
(71) (*) < exp 



-Rd(^)(0,0) 

(*) = r)/n P i N D(z) (<t>) = 0}- 



n*) . . («i) 2 1 



5^(a r 



n>l 



y(z) 3 i 

7T + 2 ~ 2^2 



tin tfo 



B(R) 2 



Prom (57), (71), notation of Theorem 4 and (60), it follows that 



E[V k {C\B{R))\ 



(72) 



(C\B(R))* 
< (3/2)-J(R) 



K(0,0) P {N D{z) (d>) = 0}dz 



/ exp 


\ V(z)l 


TT 


l(C\B(R)) k 



dz 



EV k (C p \B(R))-e 



(3/2)-J{R) 



which proves part (b) of Theorem 4. 

Problem 4. The facts that EV{C) = EV(C P ) = vr and that the Ginibre- 
Voronoi tessellation is more regular than the Poisson-Voronoi tessellation 
suggest the conjecture (which seems to be confirmed by Monte Carlo simu- 
lation [16]) that the inequality, 



EV 2 (C) <EV 2 (C 



holds. It would be interesting to provide a rigourous proof of this property. 
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5.3. Proof of Theorem 5. By Theorem 1 and formula (54), the typical 
cell C coincides in law with the cell C(0) related to the process </>o which is 
obtained by removing from <j> the point Z. If we consider the cell, 

C = {z G C; Vv G 0, \z\ < \z - v\} C C, 

that is, if we do not remove the point Z, then for every Borel set A C C and 
fc>l, 

E[V k (C nA)]= [ P{iV D(2) (0) = O}dz 

JA k 

(73) = / exp{-^l-£l / K^u)du\dz 

JA ^ n>2 J > 

< [ exp{-V(z)/Tr}dz = E[V k (C p nA)}. 
JA k 

Say that a point u G is a neighbor of the origin if the bisecting line of the 
segment [0,u] intersects the boundary of the cell Co- 

Denote by N = N(4>) the set of neighbours of the origin. Recall property 
(56), that is, 

zeC{0) k £>(z)n0 o = 0. 

By Theorem 1, we have also 

zeCn D(z)ri(t) = 

(74) 

<==>- D(z)C\<j)Q = and Z<£D(z). 

Moreover, \i Z qLN then obviously C(0) = Co- Consequently, we obtain 

E[V k (C(0))-V k (C )} 

= E[{V k (C(0))-V k (C )}xl {ZeM} ] 

= [ [P{D(z)nc/ )o = 0,Z€M} 
Jc k 

- P{D(z) r\<f)o = 0,Zi D(z), Z G M}] dz 

= I P{D{z)C\<f>Q = 0,Z G D(z),Z eAf}dz 
Jc k 

<P{ZeM} l/2 [ P{D(z)n<t> o = 0,z eD(z)} l ' 2 dz. 

Now by (53), (19), (31) and notation (5), 
P{D(z)C\4> Q = 0,Z eD(z)} 



(75) 
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= ^__K(0,0) 

(76) 



^jP{N D{z) ( ( f> ) = 0} 



V -Rd( z )(0, 

(nR D{z) (0,0)-l)P{N Diz) ^) = 0} 
= H(z). 

Moreover, by (19), 

E[V k (C(0)) - V k (C )] 

(77) =/ (P{N D{z) ^ ) = 0}-P{N D{z) (4>) = 0})dz 

H{z)dz. 
Relations (75)-(77) imply (6). 

5.4. Theorem 5, case k = \. Fix a S C and consider the kernel K(z\,Z2) 
K{z\ — a, Z2 — a); hence 

f>, ^ ^ e-^' 2 ^-^ 2 ( Zl - a) n 



n>o VTrn! 
(78) 

e -(l/2)| Z2 -a| 2 ( i — 
X ■ — 



Observe that the functions 

f n {z) = e -(V2)l*-«l a ( z -a) n , z £ D(o),n > 1, 
are orthogonal on -D(a) = B(a,r) with r = |a| and that 

|/„(z)| 2 dz = vr 7 (n + l,r 2 ), 

D(a) 

where 

-tj.n—1 



7(n,«)= / e-H n ~ l dt 
Jo 

is the incomplete gamma function. 
Denote 

7(n + 1, r 2 ) 

79 «n= rV , ' ; , n>0, 

n! 

and 

-2^-1/2. 



/ n = (7r 7 (n + l,r 2 ))- 1/ 7 ] 
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then 
(80) 



K(z 1 ,Z 2 ) = y^a n fn(zi)fn(z2)- 



n>0 



It follows from (80) that a n , n > 0, are the eigenvalues of the integral kernel 
K on D{a) = B(a,r) with r = \a\ and that for the resolvent kernel Rd(o)i 
we have 



^(a)(0,0)=^-^/ n (0)/ n (0) 

w ^ 1 - a„. 



(81) 



where 



n>0 

-y- 

7T 4^ T(n+ l,r 2 )' 

n>0 v ' ' 



r 2n e -r 2 



T(n, u) = T(n) — 7(n, u) 



e~H n - l dt. 



By Remark 17 and formulas (14), (31) and (79), 



(82) 



mw0)=o}=^ £^! ^r : 



r = a . 



n>0 



Therefore, by (81) and (82), 

W(C)= / P{iV D(2) (0o)=O}(iz 



7T 



7T 



n 



^d W (0,0) 
r(n + l,t) 



n 



x 



n>0 n>0 
t=+oo 



r(n + l,t) 



r// 



n 

n>0 



r(n + i,t) 



7T. 



t=0 



This is the expected result. Now, we have also 

EV(C )= [ P{N D{z 0) = O}dz 



7T 



+oo 



n 

n>0 



T(n + l,t) 



(It. 



IV. 
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In [20] [formula (15.1.27)], M. L. Mehta showed that, for every t > 0, 
(83) 



n £(n + l 1 t)^ (1 + t)e _ 2 , 



n>0 



7?.! 



This implies 



3-7T 

EV(C ) < —. 



From (79), (81), and with the notation (5), we obtain 
H(a) = (7rB D(a) (0, 0) - 1) J] (1 - a„; 



n 



T(n + 1, 



r 



)V. 



X 



n>0 



E 



r e 



(84) 



n>0 n>0 
..2\ 



r(n + l,r 2 ) 



n 

n>0 



r(n + 1, r 



n 



n 

n>0 



T(n + 1, 



n! 



E 

n>l 



r e 



r(n + l,r 2 



Consequently, 
(85) 



/ H{z)dz = 2Tr 1- / T ■ 
• /: ' V Jo n>0 



r(n+l,t) 



n! 



and inserting (84) and (85) in (6), we obtain (7). Now, applying the Holder 
inequality, we get 



E 



n c — t 



t n e 



n r " ,+ . M) * 



< 



r(n + l.t) I " n\ 

v n>l v ' y/ n>0 

+oo / t n e -t \ r(n+l,t) 

2^rf n + i,t)Hl 



n>l v ' ;/ n>l 



;?.! 



e~* dt = l. 



This and inequality (7) give 



P{ZeAf ((/>)}> 



n 

n>0 



r(n+l,t) 



n! 



and by (83) we obtain 
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Problem 5. It would be interesting to investigate other geometric char- 
acteristics of C, among others, the number of sides, the perimeter, and the 
radius of the smaller disc containing C (0) . 



APPENDIX: PROOF OF LEMMA 13 
Let a, P be point processes on E such that for every decreasing event 

(86) P{aeA}<P{peA}. 

Recall that C T is the collection of sets which are a finite union of ele- 
mentary decreasing events. We want to prove that point process a G ,M CT (E) 
stochastically dominates the point process f3 G 7W,j(E) which is equivalent 
to the fact that inequality (86) above is satisfied for every decreasing event 
AeT{M). 

A.l. Step I. It suffices to prove that (86) is satisfied for every decreasing 
event A G .F(E) by assuming that A G .F(E) is a closed subset of A^ CT (E). 
Moreover, if E is compact then we may assume that A G .F(E) is compact 
as well. 

Indeed, denote by Q the law of point process a and by Q' the law of 
the process /3. We want to show that Q(A) < Q'{A) for all A G T(E). Re- 
call that A^o-(E) is a Polish space. Then, by the Lusin theorem, Q(A) = 
sup{Q(2t);2l C A and 21 is compact} and Q'(A) = sup{Q'(2l); 21 C A and 21 
is compact}. Consider now a compact set 21 C A and denote 

2t = {£ G M ff (E); 3C G 21 such that f C f }. 

The set 21 is decreasing and 21 C 21 C A. Consequently the result follows from 
the lemma below. 

Lemma 26. (i) The set 21 above is closed. 

(ii) IfMis compact then the set 21 is compact as well. 

Proof. For property (i), consider a sequence (£ n )n C 21 such that ^ n — >• 
^ G M a (E) [the space M. a (K) being endowed with the vague topology]. We 
want to show that £ G 21. For every n > 1 there exists £ n G 21 such that 
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£n Q Cn- With the set 21 being compact, there exists a convergent subse- 
quence Cn k — > C £ 21. We claim that £ C £ (and thus £ € 21). Indeed, suppose 
that there exists x G £ such that x ^ £. Let :E — >• M, i = 1,2, be contin- 
uous functions with compact supports, respectively, K{, i = 1,2, such that 
#i C tf 2 , £ n = x, 1T 2 n C = 0, < h,f 2 < 1, h{x) = 1 and / 2 = 1 on 
A'l Denote F(£) = {r?; |l-E*e,/i(*)l < 1/2} and F(C) = {??; lE^WI < 
1/2}. For large we have £ nji , G V(£) and £ nfc G V(C) from which follows 
that £ nfe n K\ ^ and £ nfe n iiCi = which implies in turn the contradiction 
£,n k Cn k - Property (i) is then proved. To prove property (ii) notice (see [14]) 
that the compactness of the sets E and 21 implies that there exists A > 
such that 21 C {N E < A}. Obviously 21 is included in {iV E < A}, the later 
being compact (see [14], page 33); the result follows from property (i). □ 

A. 2. Step II. We may suppose that E is a compact separable metric 
space. Indeed, we have E = (Jn>i ^ n w ^ eie the sets K n f E are compact 
with countable bases. Denote a n = ad K n and /3 n = (3 n K n . Condition (86) 
implies that P{a n G .4} < P{f3 n G .4} is satisfied for every decreasing event 
A G J-d(K n ). Thus if condition (86) implies stochastic domination for E 
compact, then the process f3 n is stochastically dominated by the process a n , 
and we have P{a n G .4} < P{/3 n G A} for A G T(E). This and the fact that 
« n ta and j3 n "f /3 implies that P{a G A} < P{(3 G .4} for closed decreasing 
sets A£ T(E). 

A. 3. Step III. We suppose that E is compact with a metric d. Fix a com- 
pact decreasing set A G F(E). We can suppose that there exists e > such 
that for each A G A and each x,y G A, we have d(x,y) > e. Indeed, denote 
by A n C A the set where the elements are the finite sets A £ A such that 
for every x,y G A, we have d(x,y) > 1/n. Then the sets A n are decreasing 
as well and when n — > +oo, we have Q(A n ) — > Q{A) and Q'(A n ) — > Q'(A). 

A. 4. Step IV. For each A G A (note that A is finite) and n> 1/e consider 
the set 

O n , A = {0€ .Mo- (E) ; iV B ( x i/ n ) (0) < 1 for each x £ A 

(87) 

and A^\ UxeA B (x,i/n)(<£) = 0}, 

where Bq{x, 1/n) is the open ball. 
Denote Sin = LUeA °n,A- 
In order to finish it suffices to note that: 

- The sets O n> A are open (see [14]); 

- We have $in+i C Sin', 

- f)On = A; 
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- &n is a covering of A by open sets, A being compact there exists a covering 
of A by a finite number of sets O n< A- 

Consequently, in order to obtain Q(A) < Q'{A) it suffices to have 

QU. =i,...JV On,Ai) < Q '(Ui=i,...,jv On,Ai) ■ Naturally, Ui=i,...,jv O n ,Ai G which 
completes the proof. 
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